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Abstract —Two-Way Relay Channel (TWRC) plays an important role in 
relay networks, and efficient relaying protocols are particularly important 
for this model. However, existing protocols may not he able to realize 
the potential of TWRC if the two independent fading channels are not 
carefully handled. In this paper, a Hybrid DeModulate-Forward (HDMF) 
protocol is proposed to address such a problem. We first introduce the 
two basic components of HDMF - direct and differential DMF, and then 
propose the key decision criterion for HDMF based on the corresponding 
log-likelihood ratios. We further enhance the protocol so that it can be 
applied independently from the modulation schemes. Through extensive 
mathematical analysis, theoretical performance of the proposed protocol 
is investigated. By comparing with existing protocols, the proposed HDMF 
has lower error rate. A novel scheduling scheme for the proposed protocol 
is introduced, which has lower length than the benchmark method. The 
results also reveal the protocol’s potential to improve spectrum efficiency 
of relay channels with unbalanced bilateral traffic. 

Index Terms —Two-way, relay, hybrid demoduiate-forward, maximum 
likelihood detection. Log-likelihood ratio, queue analysis, scheduling 

I. Introduction 

A typical Two-Way Relay Channel (TWRC) involves two source 
nodes (A and B) exchanging data through one intermediate relay 
(Relay) m-oi, as shown in Fig[T] The two-way problem was 
first studied by Shannon {4) and analyzed theoretically in terms of 
capacity by O. Recent research confirms that the employment of 
relays in certain conditions can significantly increase the performance 
of wireless networks if they are handled effectively using suitable 
relaying protocols 0, (3 To this end, several common relaying 
protocols, e.g., Amplify-Forward (AF) and Decode-Forward (DF) 
were proposed. However, AF has the problem of noise amplification 
and involves expensive RF chains to mitigate the existing coupling 
effects, thus practical DF protocols have been proposed which, 
through proper demodulator design, can also have better performance 
and achieve full diversity (8). 

The TWRC model has extensive real-world instances, such as 
ground stations connected by satellites, two mobile users exchanging 
information through one common base station or relay, or sensors 
connected by a smart router, for example ©-ini- In our recent 
project on developing UAVs to bridge communications between 
wireless users, such a model can be directly adopted 03, HU- 
As the two source nodes can transmit signals simultaneously, an 
effective relaying protocol with high efficiency is critical to system 
capacity. However, the application of some simple protocols, e.g. 
AF/DF, incurs low bandwidth efficiency because four time slots are 
required to finish the exchange of one message from each source. 

Considering the special structure of TWRC, researchers have 
applied the concept of network coding in the wireless domain and 
proposed several protocols that need only two or three time slots. 
As a result, significant spectrum usage (or time)-at least 50%-can be 
saved. Such protocols include Analog Network Coding (ANC) CD, 
Physical-layer Network Coding (PNC) GD. Digital Network Coding 
(DNC) (H and their equivalences. Even though these new protocols 


can theoretically improve system throughput by up to 100% over 
the direct relaying transmission scheme, due to the limit of channel 
effects such as fading, Signal-to-Noise Ratio (SNR) variation and 
system implementation difficulties such as imperfect synchronization, 
the claimed advantages cannot be fully achieved in real-life scenarios. 
For example, if channel gains of the two source-relay channels have 
large discrepancies, these schemes may incur large packet error rates 
and spectrum efficiency is greatly decreased by retransmitting large 
amounts of erroneous packets Qsi, mi- This problem is particularly 
severe for the PNC protocol as it achieves the best performance when 
both channel gains are the same. It is therefore of theoretical interest 
and practical importance to design protocols/network coding schemes 
which can overcome the problems of wireless channels with dynamic 
fading and large SNR variations. This is the problem which this paper 
wishes to solve. 

Current efforts on the differential demodulation of received super¬ 
posed signals to reduce noises have been substantially studied in the 
literatures Qa-m]. Such modulation is an example of the compute- 
forward protocol used for PNC. The authors of m summarized the 
necessary and sufficient conditions on the mapping of symbols at 
the user side to remove the ambiguity at relay, and thus improved 
the error rate performance of a system adopting phase-shift keying 
modulations. A novel DeNoise-and-Forward (DNF) scheme was 
proposed in (m, which addressed the noise accumulation issue 
at the relay node of a two-way system and inspired the proposed 
modulation-independence protocol of this paper. The analysis from 
1201 shows that at the source to relay phase, the traditional operation 
of PNC - ExclusiveOR does not always offer the best performance, 
thus the authors proposed a design strategy to optimize constellations 
for improved throughput. Similar focus on the constellation design 
can also be found in HD which aimed at removing ambiguity points 
of M-PAM (pulse amplitude modulation) signal constellations to 
enable using binary physical-layer network coding in an efficient way. 
This paper also provided an instrument for higher order modulation 
schemes. Even though extensive studies have been carried out, 
practical scenarios with variable channels and large SNR variations 
with the focus on low error rate and queue scheduling are still yet to 
be explored and leave the gap for further work. 

In this paper, we propose a Hybrid DeModulate-Eoward (HDMF) 
protocol specialized for the TWRC model. Specifically, the HDMF 
protocol focuses on two-way relay channels which do not always have 
equal channel gains because of fading and synchronization errors. 
The aim of this protocol is to achieve both high spectrum efficiency 
and low packet error rate when compared with the typical protocols 
introduced previously. HDMF incorporates differential and direct 
DMF, and adapts its strategy based on the channel coefficients to 
reduce packet errors. A key to this protocol is the decision rule which 
regulates the two DMF schemes. The corresponding Log-Likelihood 
Ratios (LLR) are used to form the decision criterion. We also improve 
the simple HDMF protocol for the case where high dimensional 
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Fig. 1. 

Two-way relay channel model. 



modulation schemes are used. The proposed HDMF does not only 
increase spectrum efficiency over direct relay schemes, but also 
overcome the limits of ANC/PNC schemes. The proposed HDMF 
can be directly extended to the hybrid-DF case, however, since this 
only increases the complexity (by adding a decoder) without bringing 
more understanding about the system, hybrid-DF is not covered in 
this paper. 

Furthermore, we propose a scheduling scheme to adjust the data 
flows of HDMF considering dynamic channel conditions. We in¬ 
troduce four transmission modes and develop a scheduling scheme 
based on the transition probabilities of the modes. Without loss of 
generality, we only study the transition states of the queue length at 
source node A and the forwarding queue length from A to B at Relay. 
Similar behaviour can be observed at B and also the queue at Relay 
(from B to A), thus they are omitted in the paper. Queue behaviour 
is summarized for all four modes and is an essential component of 
the scheduling scheme: The average length of the queue at A and 
the forward queue at Relay for data to B are investigated based on 
a Markov chain model. Its performance is verified using simulation 
and compared with the DNC scheme studied in Gl. 

The main contributions of this paper are summarized as follows. 

• A HDMF protocol proposal for the TWRC model with fading 
channels and variable SNR; 

• A scheduling scheme proposal for the TWRC model with the 
HDMF protocol; 

• A study of average queue lengths of the scheduling scheme with 
the HDMF protocol. 

This paper continues as follows: Section II describes the system 
model; Section nni proposes the hybrid DMF protocol; Section lYl 
designs the scheduling scheme and investigates the queue length 
transition states; Section lYll verifies its performance by simulation; 
Section ivni discusses the important issues; Section IVIIII concludes 
the paper. 


H. The System Model 

The system model (Fig. [TJ considered in this work is the typi¬ 
cal TWRC model with two source nodes to exchange information 
through an intermediate relay. The channels are assumed to remain 
unchanged for at least one packet duration and the coefficients are 
known by the corresponding destinations. In designing the protocol, 
we assume noise to have a Gaussian distribution with zero mean and 
N(,/2 variance. 

At every time slot, sources A and B simultaneously transmit 
data packets, e.g., Xa(n) and Xi,(n) to Relay. Each packet has 
M symbols, denoted as Xa(n) = [Xa{n),--- xt,{n) = 


[xl(n), ■ ■ • , x ^Relay receives superposed versions of the two 
packets from A and B as follows 

Yrin) = /tAJt(n)xa(n) + hBR{n)xb{n) + Wr(n), (1) 

where yr(n) = [yj.(n), ■ ■ ■ , (n)]'^, hAR(n) and hBR{n) are the 

channel coefficients between A and R, B and R at the time slot n. 
Wr(n) is the noise vector at Relay. The mth symbol in yr{n) can 
be denoted as y^(n). 

Messages received at Relay are processed and forwarded to the 
corresponding destinations by protocols, such as DMF, AF and DR 
The forwarded signal Xr(n) is given by 

Xr(n) =/(y^(n - 1)), (2) 

where yr{n — 1) is the received signal by Relay at the (n— l)th time 
slot and /(■) is the relaying function which describes the processes 
of the relaying protocol. For example, if using direct DMF 1171 . /(•) 
denotes the process of demodulating the received signal and then re¬ 
modulating it using the designated modulation scheme. In this case, 
data from the other user will not be detected. However, such case 
only happens if the channel quality deteriorates to a considerably low 
level as denoted by the likelihood ratio value. In order to ensure the 
dropped data can be retransmitted, current communication systems 
usually adopt Automatic Repeat reQuest (ARQ) schemes which are 
also assumed in this paper. The implementation details of ARQ are 
omitted here and can be found in dn. If the weaker channel becomes 
strong enough, differential DMF will be automatically switched, 
where /(•) denotes the demodulation and remodulation of the data 
from both users. Even though sequential interference cancellation 
techniques can be used, it is not the focus here since the error rate 
of the weaker channel would be even higher than the stronger one. 

The received signals at source B and A at the nth time slot are 
given by 

yb{n) = /ins(n)xr(n) -f Wb{n), 
ya(n) = ftnA(n)xr(n) -I- Wa(n), 

where yb{n) = [ylin),--- hRB{n) is the channel 

coefficient between R and B at the nth time slot, and W6(n) is the 
noise vector at node B. ya{n) and hRA{n) have similar definitions 
as those for source B. 

HI. Hybrid Demodulation-Forward Protocol 

This section proposes the HDMF protocol. The general idea of 
HDMF is to have two DMF modes including direct DMF and 
differential DMF, and switch the two modes automatically based on 
a criterion obtained from the channel gains and SNRs. The protocol 
is implemented on the Relay and destinations respectively, where the 
Relay uses it to construct and forward signals. After the reception of 
these signals, A and B apply the protocol to detect the source data. 
We use Quadrature Phase Shift Keying (QPSK) as an example to 
explain the mechanism of this protocol which can be easily extended 
to higher-order modulation scenarios. The details are elaborated in 
the final part of this section. 

To introduce the principles of this protocol: we assume that data 
packets from the two source nodes are of the same length. (The 
case of unequal length packets will be discussed in Section IWl ) 
Each packet has M symbols, and each symbol denotes K bits if the 
modulation scheme has K-th order. For example, the source node 
packet Xa(n) has M symbols and the mth symbol x^(n) has K bits 
which are denoted as {fe2^i(n), • • • Similarly, x^{n) has 

K bits and are denoted as {6^i(n), • • • We assume the 

symbols are generated randomly with equal probabilities. Therefore, 
for Quadrature Phase Shift Key (QPSK) modulation, each symbol has 
two bits since K = 2, and the four different constellation symbols 
of QPSK are generated randomly with the same probability: i. 
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Fig. 2. The demodulation scheme of direct DMF where 
and QPSK. The dark shadow areas around the two decision lines denote high 
detection error prohahilities. The corresponding symbols and data bits are 
mapped using Gray code and shown in the four corners. 


The following subsections first introduce the two component parts 
of the HDMF protocol, and then present the decision criterion to 
choose these two components and the implementation details. 

A. Direct Demodulation-Forward 

The received signal at Relay is shown in Q. As two signals arrive 
at Relay simultaneously, we should try to demodulate the stronger 
signal directly in order to control the detection errors if one channel 
has significantly higher channel gain than the other and the noise is 
strong. Without loss of generality, we take the route from A to Relay 
to B as an example. 

The maximum likelihood (ML) criterion is used to demodulate the 
mth symbol within Xa(n), 

x^{n) = a.rg majc {P{y^{n)\x'^{ri)} , rri = I,M. ( 4 ) 

where AA is the modulation symbol set. For QPSK, Af = {1 + 
i, —1 + i, —1 — i, 1 — i}. Therefore, direct DMF has the forwarding 
signal as follows 

Xr(n) = Xa(n) = ^xl{n), xf^(n) j (5) 

In order to describe the scheme concisely, the symbol index m is 
omitted henceforward and Xa always denotes the mth symbol within 

Xa. 

Since all symbols are drawn under the same probability, from Q 
we have 

P(yr-(n)\xa{n)) 

= ^ P{xb{n)) ■ P{yrin)\xain),Xb{n)). ( 6 ) 

Given QPSK modulation, P{xb{n)) = 1/4, Xbin) £ M. ® can 
be simplified as 

P{Vr(n)\Xa{n)) 

= i P{yr{n)\xa{n),Xb{n)). 

Under the ML criterion, the symbol (Xa(ri)) with the maximum 
P(yr{n)\xa{n)) value is chosen to be the detected symbol within 
Xr(n). 
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Fig. 3. The demodulation scheme of direct DMF under BPSK modulation. 


Fig. [2] shows the decision areas of the direct DMF scheme for the 
case of |/iAfl| > \hBR\ under QPSK modulation. The other case 
when |/iah| < |hsij| can be similarly analyzed and is omitted here. 
In Fig|2 the decision lines are the axes. The detected symbols and 
the corresponding data bits are given in the tables. The areas are 
derived based on the distances between the received signal position 
and the decision lines for illustration purposes. (The decision areas 
for Binary Phase-Shift Keying (BPSK) modulation is also given here 
and shown in Fig[3 which is a simpler version of the QPSK case.) 

Ideally, the received symbols should be located at the four con¬ 
stellation points in the quadrants of the Cartesian system. However, 
because of noise and interference from the other channel, the received 
signals of y^. are scattered on the plane around the four lines 
(|hAfi(n)| ± \hBR{n)\, |/iAfl(n) ± hBR(n)\ x i) which form the 
two stripes in the figure. If some received symbols fall into the dark 
areas, a small noise can push them to the other side of the decision 
line, which will then be interpreted as different symbols and cause 
detection errors. For symbols located in the lighter areas, because 
they are far away from the decision lines, they need stronger noise 
(or lower SNR) to cause demodulation errors. 

Since the received signals are around the above four lines, in order 
to have less detection error, the lines should be far away from the 
axes so that even though some symbols are located in the dark areas, 
they would still be far from the decision lines and are strong enough 
to combat noise. Therefore the width of the dark stripes is essential 
for direct DMF. From Q, we can see that this width is controlled 
by the relevant channel gains of the two channels - the larger the 
discrepancies between the two channels, the wider the stripes in the 
figure. The ideal case to apply direct DMF is when one channel is 
much stronger than the other so that we have large absolute value of 
\hAR{n)\ — \hBR{n)\, thus the lines are far away from the axes. 

By way of contrast, we can see where direct DMF will introduce 
problems by looking at an opposite example where the two channels 
have equal channel coefficients. The stripe width will be 0 and the 
received symbols will be around the axes (which are the decision 
lines) and the lines with the value of double channel gains. We can 
therefore expect | of error rate because a small noise can push 
the received signals around the axes from one side to the other 
side. The clear boundaries where direct DMF should be applied are 
essential to the HDMF protocol and should be jointly considered with 
SNR. Before we discuss this subject, we introduce a solution to the 
scenarios where the direct DMF fails. 


B. Differential Demodulation-Forward 

Previous discussion reveals that direct DMF would generate results 
with high error rates if the two channels of TWRC have comparable 
gains, e.g. \hAR{n)\ « |hfli?(w)|. In order to solve this problem, 
we propose to use differential DMF rather than direct DMF in some 
scenarios. Fundamentals of this scheme are introduced using QPSK 
as an example. 

1) Differential DMF at Relay: In order to demonstrate the differ¬ 
ential DMF algorithm, we firstly consider an instance where there is 
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no noise and hARin) = hBR{n). In this case, the two source signals 
of Q arrive with the same power at Relay. The received superposed 
signal can be denoted as 

yr{n) = hAR{n){yia{n) + X6(n)). (8) 

For a given modulation scheme like QPSK, both the real and 
imaginary part of yriji) only have two possible absolute values: 
2|/iAfl(n)| and 0 respectively. These two values correspond to the 
two input sets: ba,kin) = h^kin) and ba,kin) / 65 ,*(n). It is 
easy for the relay to directly demodulate the received symbols at 
the first instance as both of them are the same and the power is 
enhanced. However, the latter case adds difficulty as it may be the 
superposition of two different bit sets: ha,kin) = 1, bh,k{n) = 0 
and ba,k{n) — 0, bf,,k{n) = 1, (fc = 1,2). However, since both the 
two sources know their transmitted signals, it is desirable for both A 
and B that Relay differentially demodulates and forwards signals to 
them. 

The differential DMF at Relay, after converting data symbols to 
data bits, is given as 

br,k{n) = ha(sb,k{n) = ba,k{n) © bb,k{n), fc = 1, 2, (9) 

where © is the XOR (exclusive OR) operation. For example, if node 
A has a symbol 1+i and node B has -1+i, the corresponding bits of 
A’s symbol are 1,1 and those of B’s symbol are 0,1. The resulting two 
bits after the above processing can be given as br,i{n) = 1 © 0 = 1 
and br, 2 {n) = 1 © 1 = 0 . 

The above signal processing can be applied directly in the symbol 
domain if we can demodulate the differential between ba,k{n) and 
bb,k{n) as follows 

fca06,A;(rt) = 

f 0, Re{i/r(n)} = 2Re{/iAfl(n)}; 1, Re{yr(n)} = 0 (10) 

\ 0,lm{yr(n)} = 2Im{fcAfl(n)}; 1,Im{i/r(n)} = 0 

where Re{-} and Im{-} denote the real and imaginary part of the 
input, respectively. The first line of (Uni is for the first bit and the 
second line is for the second bit because one QPSK symbol denotes 
two bits. 

The ML differential detection at relay is given by 

XaSibin) = a.Tg max {P{yr{n)\xa(sbin)}, (11) 

where P{yr{n)\xa^b{n) G Ad) is the distributions of yr{n) given 
the differential symbols within modulation set M. 

From m, it can be seen that 

Piyrin)\xafsbin)) = ^ 

Xa('n-)€J<4,xij(n)GA4 (12) 

P{xa{n),Xb{n))Piyrin)\xain) © Xbin)). 

Under the conditions of QPSK modulation and symbols generated 
with equal probability, 

Piyr[n)\Xa(Bb[n)) 

= \ Piyrin)\xa{n)®Xbin)). 

Xa('n)GJ^,Xf,(n}GAi 

Relay can demodulate the received signals differentially using the 
above ML detection algorithm and forward the results to A and B. 
Take the case of |fcAJt| > Ifcsjtl for example. The demodulation 
scheme of differential DMF is shown in Fig. [4] Similar to Fig. 
the dark stripes are the areas which has more detection errors than 
the rest. The decision lines are formed by ©I/iahI and ±|/iAfl|*. The 
edges are formed by ±|/iAfl|±|fcBfl| and their imaginary equivalents. 

Similarly, the received symbols are located around the lines which 
form the four stripes. However, in contrast to direct DMF, the decision 
lines are not the two axes anymore. Instead, they are ±|/iah| and 
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Fig. 4. The demodulation scheme of differential DMF under QPSK 
modulation. If the received symbols are located at the dark shadow areas, 
they have potentially high detection errors. (Assume ^ \^br\)- 
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Fig. 5. The demodulation scheme of differential DMF under BPSK 
modulation. 


©Ifc-Afll*. In this case, the symbols which are located around the 
axes can mostly combat higher noise. The ideal case is when the 
two channel gains are the same, e.g. {HarI = Ifc-BJil, (where direct 
DMF has the worst performance) so that the majority of the received 
symbols will be located around the axes (not depicted in the figure) 
and far away from the decision lines. In this case, they will strongly 
combat noise. The width of the stripes will decrease if the two 
channels have different gains. The worst case is when one of them is 
zero, e.g. |/isb| = 0 in Fig. |4] The width, in this case, will become 
zero and the edges (mean of the received symbols) overlap with the 
decision lines. In this case, a small noise will push the signals from 
the correct side to the wrong side, thus generating a large error rate 
of |. Another important job is to decide when to use differential 
DMF. (Fig|5]gives the decision areas for the similar case under BPSK 
modulation, from which we can see when yR has values around the 
two decision lines: — |/iab| and |/iab|, the detected symbol has a 
higher error probability.) 

2 ) Detection of differential DMF symbols at A and B: The signal 
received by the two terminal nodes A and B are given by (O. We 
use node A, for example: ya(n) = fcijA(n)xr(n) + Wa(n) and 
Xr(n) = f(yr(n — 1)) (given by ([2l(). At the nth time slot, if the 
differential DMF protocol is applied to generate Xr{n), we obtain 
br,k{n) = ba^b,kin— 1) for the fcth bit of the mth symbol of Xr(n). 
These received symbols should be demodulated first and then detected 
by the differential DMF protocol. 

Similar to ML detection used in direct DMF, the differential DMF 
symbols received at the destination are firstly detected using (|4j. After 
the ML detection, we have the estimated symbol Xr(n). The equation 
of this step is similar to (0 and is neglected here. 

The second step is to decode the original data bits of the source 
symbols Xain — 1) and Xb{n — 1) through an XOR operation as 







































5 


follows 

ha,k{n) = hh,k{n - 1) © hr,k{n), fc = 1, 2, 

ib,k{n) = ba,k{n - 1) © br,k{n), fc = 1,2, 

where br,k{n) is the fcth data bit of the estimated symbol Xr{n). The 
other parameters of the equations are ba,k{n — 1) and bb,k{n — 1), 
which are the source data bits transmitted by the corresponding nodes, 
respectively, at one time slot previously. Such data should be stored 
in source nodes for later use. We can then obtain the desired data 
through the XOR operation in ( 114b . 

Comparing the error areas of the two DMF schemes, we can 
find that differential DMF and direct DMF form a complementary 
relationship: one excels when the other has poor performance. It 
is easy to know which one should be used for the extreme cases 
discussed before. However, for a practical fading TWRC model 
whose channels are fading, the scenario is usually between these 
extreme situations. Therefore, it is important to find the decision 
criterion for the appropriate application of HDMF. 


communication quality. Thus the decision rule for a packet to be 
processed through direct or differential DMF is given as follows 

mm{\Cuir{n)\,m = 1,..., M} > mln{|£Si/(«■)!, 
m = l, ...,M}, Direct DMF, 
min{|£Sir(tt)|,tTi = 1,..., M} < min{|Z:Si/(n)|, (19) 

m = 1,..., M}, Differential DMF, 

\C^Un)\ = max{|£St.A(n)|,|£S,.s(n)|}. 

We select the smallest LLR value from one packet transmitted from 
the stronger channel under direct DMF and compare it with the 
smallest LLR under differential DMF. By using this criterion, the 
DMF scheme with a greater minimum value will be selected to 
forward data. Based on PER, even if there is only one failing symbol, 
the whole packet is labelled as in error. The criterion above therefore 
ensures that the one with lower PER is selected. 


D. HDMF Detection at End Nodes 


C. Implement HDMF at Relay 


As discussed above, direct and differential DMEs are more suitable 
for different channel settings of the TWRC model. A key problem is 
to apply the ML decision criterion to choose the right DME schemes. 
The Log-Likelihood Ratio (LLR) tool can be used to implement such 
criterion. 

The LLR value of each bit of the symbols under direct DMF 
is calculated based on channel conditions and SNR of the received 
signals, as follows, 


CoirAin, k) = log 


C.DirB{n, k) = log 


P(yr(w)|ba,fc(n) 

P{yr{n)\ba,k{n) 

P{yr{n)\bb,k{n) 

P{yr{n)\bb,k{n) 


JO 

0 )’ 

JO 

0 )’ 


(15) 


where jCr>irA{n, k) is the LLR of the fcth bit of the mth symbol 
from A at the current time slot and jCoirB{n, k) is that obtained 
from B, where the index m is omitted in the equations for simplicity. 
Similarly, the LLR value under differential DMF can be obtained as 
follows. 


Coifin, k) 


, P(yr{n)\ha^k{n) 7 ^ bb,fc(n)) 
P{yT{n)\ba,k{n) = bb,k{n)) ’ 


(16) 


where €.oif{n,k) is the LLR of the fcth bit under the protocol of 
differential DMF given Xain) and Xb{n). 

It is easy to know that the sign of an LLR value under the 
ML criterion denotes the detected bit result and its absolute value 
denotes the degree of confidence (23). For a packet containing several 
symbols, we can try to find the symbol with the minimum confidence, 
which would most likely cause high Packet Error Rate (PER). 

The LLR values of all symbols (where each symbol has K bits) 
in the received packet at Relay, e.g., yr{n), (M total symbols) under 
direct DME are calculated as follows. 


P'DirAijt) — ^ ^ \F'DirA{cii fc) 1 7 


K 


(17) 


PDir bIcI^ — ^ ^ \I~-Dir b(xL^ k^\. 


Similarly, the LLR values of differential DME are calculated as 


K 

C-Difin) = ^ |£oi/(n, fc)|. 

fe = l 


(18) 


The Relay with HDMF protocol can generate and forward two 
different kinds of messages by using the direct or differential DMF. 
For the end nodes (A and B), it is essential to be able to detect which 
kind of messages are transmitted from Relay. A possible solution is 
to add one bit in the packet header at the Relay, which indicates the 
DMF scheme used. However, such method increases implementation 
complexity and changes the frame structure, which is less desirable 
for practical wireless systems. In this section, we propose a blind 
detection algorithm with low complexity. 

At the end nodes, e.g., source B, the received packet yb{n) is 
first demodulated and mapped to data bits. A Cyclic Redundancy 
Check (CRC) is then performed on these data bits. If it is correct, 
this message is processed by direct DMF at Relay. If it is not correct, 
the second detection attempt using differential DMF (See Section 
IIII-B2b will be carried out. The decoded data bits, e.g., ba,k{n), are 
checked by the CRC. If it is correct, differential DMF is used in the 
relay. Otherwise, the packet has some error symbols and should be 
discarded. 

Finally, we analyze the above detection algorithm to see whether 
it can function properly in real life conditions. For a packet with a 
large enough number of symbols, e.g., M = 128, it contains 256 bits 
under QPSK modulation. (Higher order modulation has even more 
data bits.) We assume the data bits are randomly generated with equal 
probabilities for 0 and 1. A bit level XOR operation with another 
randomly generated data packet would change approximately 50% 
of the data bits. Therefore, for a packet with 256 bits, approximately 
128 bits would be changed after the XOR operation. This would 
identify the difference between direct DMF and differential DMF 
schemes and the probability of making a mistake is quite low. (If we 
consider a simple system without source coding and channel coding, 
the error rate is approximately 1/2^^®. With coding, the error rate 
can be more than this value because of coding dependency.) 


E. A Protocol Independent from Modulation Schemes 

The hybrid DMF scheme is introduced using QPSK as an example. 
In order to apply HDMF in a system with any modulation scheme, 
the following technique is proposed: At the Relay, instead of the 
usual demodulation schemes, we introduce bitwise detection using 
the LLR values calculated in ( 115b and ( 116b . Detection of the fcth data 
bit of the mth symbol is given below for direct DMF, 


A typical digital communication system usually uses packets to 
transmit data and PER can be used as the performance index to denote 


br,k{n) = 


1, If CDir{n, k) > 0; 
0, If Coir in, k) < 0. 


( 20 ) 
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If differential DMF is selected through the criterion introduced in 
Section UlI-CI the following detection is used 


i>r,k (tl) 


1, If Coifin, k) > 0; 
0, If jCoifin, k) < 0. 


( 21 ) 


The Relay will then modulate these detected data bits into symbols 
following the designated modulation schemes and forward the results 
to the end nodes. The detection of these symbols at A and B is similar 
to the case of QPSK (e.g. using the ML algorithm) and the details 
are neglected here. 


IV. Mathematical Analysis 


This section analyzes the proposed scheme and try to establish 
the mathematical expression for the end-to-end instantaneous Symbol 
Error Rate (SER) and obtain the ratios for packets sending through 
direct and differential DMF. It is easy to know that errors usually 
come from two phases: the multiple access phase and the broadcast 
phase. In order to simplify the analysis, we take BPSK as an example. 
The error rate at relay is denoted by Pr and the error rates at the 
destination node A and B are denoted by Pra and Prb respectively. 
The end-to-end SER from A to B can be expressed as follows 

Pab^l- {I- Pr){l- Prb)- PrPrb ( 22 ) 

Similarly Pba can be obtained. The instantaneous SER can be 
expressed as 

PhDMF = - (Pab + Pba) ^ 2 ) 

= Pr P (0-5 — Pr){Pra + Prb) 


The average SER of HDMF can be obtained as 
^{Phdmf} — Phdmf 

fia)fiP)fi-y)f{0 da d^ d 7 dC 


(24) 


where a = \hAR\, P = \hBR\, 7 = IHra] and C = which 

follow Rayleigh distributions with the probability density function 
as f{x) = \ where 5^ is the Rayleigh distribution 

parameter. We can assume these channel variables are independent 
and identically distributed (i.i.d.), so that the above equation is then 
simplified to 


E{PhDMf} = E{P,} -F (0.5 - E{Pr}){E{Pra} + E{Prb}), 

(25) 

where E{Pra}, E{Prb} and E{Pr} are the average SER of the Relay- 
A channel, Relay-B channel and sources-relay channels, respectively. 

It is easy to obtain the average error rate for the relay to source 
channels as follows. 


E{Pra} = / Praf{'y)d'y 

) 

Q 




/(7)d7 


(26) 




Pra “F 1 


where pra = ^rajo^a, denotes the signal to noise ratio (SNR) and 
Q{x} = JJ” e“‘ ^^&t. Similarly 


=2(1- 


prb 


Prb H“ 1 


(27) 


The average SER of the sources-relay channels composes of two 
components from the direct demodulation (A-Relay and B-Relay) and 


differential demodulation (A&B-Relay), 

E{Pr} = [ [ {PabrPabr + ParPar + PbrPbr) /(o)/(/3)dad^, 

Ja J a 

(28) 

where Par and Pbr are the instantaneous SERs of the A-Relay and 
B-Relay channel, respectively, par and pbr are the probabilities to 
choose direct DMF. Pabr is the instantaneous SER of the differential 
DME and pabr is the probability to choose such scheme. 

Eor simplicity, l l28t is divided into two parts: the differential DME 
part and the direct DME part. 


E{Pr} = E{PDi/} -F E{PDir}, 


where 


E{PDif} = [ [ PabrPabr f {a) f{f3)dadp, 

J ct J 8 


E{PDir} = 


=// 


(29) 


(30) 


(ParPar + PbrPbr) /(a)/(^)dad^. 


Eollowing the same ML criterion used in detecting symbols, the 
probability to select differential DME ipabr) or direct DMF ( par, 
Pbr) can be obtained. Firstly, we consider the ML detection of a 
differential symbol, 

Xa(Bbin) = a.Tg max {P{yr{n)\xaq)bin)} , (31) 

the error rate of differential detection is given as 

Pabr = = l|a:r © = 0) + P(x(^ = 0|x(r © a;™ = 1)) 

(32) 

The above equation can be bounded by [ mi, eq. (7)] as follows, 

Q ( V2 min {a^Pa, P'^pb}') < Pabr 

' _ ^ _ (33) 

< Q {^^2q^Pc^ + Q (^y/2l3^pb'^ 

The same rule of GU is used to analyze the average SER as 
follows, 

pabr ~ P{<t>abr > mSLx{(j)a,4lb}} 

Par ~ ^ 4^abr © (j^a ^ 4^b} (34) 

Pbr ~ ^ (j^abr (7 0b ^ 0a} 

Since the ML criterion is used to obtain l l33t . we can derive the pa¬ 
rameters following similar mles. Given the exponentially decreasing 
characteristic of Q{x) function, the two bounds of l l33t become tight 
when SNRs are large. The lower bound of Pabr was selected in the 
calculation for its function to highlight the low-bound of the proposed 
protocol’s SER, 

Pabr ~ Q ^ ( 35 ^ 

from which we can derive the receiving SNR of an equivalent channel 
as 0a6r = 2 min{a^pa, 0^P6}- The receiving SNR of the A-Relay 
channel is 0 a = a^pa for a channel gain of c?, and that of the 
B-Relay channel is 06 = 0^P6- 

Therefore the following approximation can be obtained, 

Pabr ~ P{P^ Pb > O? Pa > ^f^^Pb © oi^Pa > /3^P6 > pa} 

Par ~ P{Q^ pa > 2)3^ pb} 

Pbr ~ p{(3^pb > 2a^pa} 

(36) 

They are used for the calculation of average SERs. 

It is easy to see that the two parts of pabr following a symmetrical 
pattern. Since all the channels are i.i.d. and the transmitting power is 
the same in all nodes, we can apply the symmetry feature to simplify 
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the calculation as follows 


POO / \ 

HPDif}~j J ^Q[V2^jf(a)f{p)dad^ 

^ V 2 pa 

+ f°° Q {V2^) mf{a)dpda (37) 

“\/ 2^Pb 

^OO / \ 

=2 y J Q /(a)/(/3)dad/3 

By introducing the full expression of Q function, we have 

^ r 


j <t<oo. Reordering these two integrals we have 
E{PDir} ~ 


POO POO 

'■i /, 


/lii v/lT" 

y' Pa 
POO POO 

"Vo 

= E{Poir} + E{Poir} 


rPt 

JbJ 


a , 

— e 272 (Jq, 

2Ph (52 
Pa 


.nil 


.izil 


‘Pa — e 2S2 


/(/3)dM/3 


/(/3)dM/3, 


(44) 


E{Pz,./} ^2 r / ^ ^ /■ 

Jo JbJ- 0 ^ Lv^Jv 


e 2 (ji 


/2paCt 


(38) 


X /(a)/(/3)dad^ 

Reordering the t wo integrals of t and a, with the new regions of 
integrals for /3^^- < a < /3y^^ snd P^/pb < f < oo, we have 

rpir 

/(a)da 


where the regions of a and t are transformed to V ^ < f < oo 
and /3^~' < Q < Pi- 

Reordering the two integrals within the above equation, the first 
term of (Ell can be written as 

E{PDip} 


POO POO 1 f 

E{Pn.f}^2 / —/ 

^0 J V^TT J j3 

POO POO ^ 

do Ji3^V2 


BVn 

=E{P'oif} -E{P'iif}. 


f Ph 

2pa 


PbP 


_ 

4(5^ Pa _ P 2(5^ p 


/(/3)dM/3 
f{P)dtdp 
(39) 


/■^ 1 (2 r poo a _^P2 p6 N 

JPlir [Jq 52 

L 


dp 


dt 


(45) 


2£Ji 2£!> + 1 


2pa 


The first term of the above equation can be calculated as 


E{P( 


u.f} =2 r f 

Jo JB 


1 -Pl P 

-e 2 J_e 2;5^V2pa^ ) 


/3v7T 


dtdp 


-Q ( 

2pb p Pa \\ Pa 
where we have used the change of variable in the first step. 

The second term of El can be similarly calculated by reordering 
the two integrals of /3 and t, 

EiPoir} 

POO 1 or POO n ^ 



= - 2 


[Jo 

1 _p 

/ __e ‘ 

J./2P6 -v/lT" 


e 25-( 


M/3 


df 


(46) 


f2 + 1 


In order to obtain a closed form of the integration, lower edge of the 
integration region is applied to obtain the following bound. 


2pa + p6 , . , I -H -—_ -I-_ 

2PaS Pb^ 

Similarly, the second term can be calculated as 

pa I . 1 


E{Pcip} > 


-2 


1 + 




(47) 


E{PBi/} = 


Pa + Pb 


1 - 


1 H-H- ^ 

PaS‘‘ Pb^ 


(41) 


For the direct DMF part, with the same assumptions of i.i.d. fading 
and the symmetrical regions of par and pbr, the two channels should 
have the same SER as follows 


E^Poir} ~ 2 / ParPardot. 


L 


(42) 


where the instantaneous SER of A-Relay channel under the ML 
criterion can be easily obtained as Par = Q{‘\J'pT ^i/p ) 1251 . 
Following the previous assumption of high SNR, an approximation 
is obtained as Par ~ Qi^) ^nd then the average SER is given as 


It is easy to see that the decreasing rate of the exponential function 
e 2 is much faster than with the increase of t. The value 

above thus provides a close approximation. 

With the availability of l EOt . (ED> ESl and ( HTl l. the average SER 
of the source to relay channel l l29t can be directly obtained. The 
results are given in Figll II 

The closed forms for the average value of the selection probabilities 
l l36t and its preceding equation l l34t provides a general measure about 
the ratio of symbols through differential DMF and direct DMF. They 
can reveal the underlying properties of the protocol. From the analysis 
in (|36) and the PDFs of a and /3, the average value of pabr can be 
obtained as 

flE 


follows 


E{Pmr} « 2y^°“ p Q f{a)f{p)dadp 


POO rP\ —— poo POt.j^-^ 

E{pabr} = / ^ f{a)f{P)dadP + / /(/3)/(a)d/3dc 

■fo J pJPi- Jo J aj-£fa 

2pb 


/o 

2pa 


2pb 


Pa 


+ 


Pb 


2pa P Pb pa P pb 2pb p Pa pa P Pb 


= 2 


P 00 PO 

Jo J \ 


POO 

1 f°° f2 

y Pa ^ 

-y= / e--df 


(48) 


/(a)/(/3)dad^ 
(43) 

The regions of integration for a and t are < cr < 00 and 


Similarly, we can obtain the closed forms for the average value of 
A-Relay channel and B-Relay channel: 


POO POO 

E{Par}=/ / _/(a)/(/3)dQd^ : 

Jo J/ 3 ,/ 2 £i. 


pa P 2pb 


(49) 




















































E{p6r} 



/(^)/(«)d/3da 


Pb 

‘2pa + pb 


(50) 


V. Transmission Scheduling over Fading TWRCs and 
Queue Length Analysis 

The two wireless channels of TWRC normally suffer from random 
fading which deteriorates the communication quality and reduces 
channel capacity if it is not well handled. Due to the random 
feature of fading, different channels of TWRC observe different 
instant fading, and it is therefore usually the case that they have 
different strengths. Under such conditions, the traditional scheduling 
schemes, e.g., the scheduling scheme proposed for DNC ca, which 
assume equal channel gains, incur large queues at Relay and decrease 
bandwidth efficiency. On the other hand, with the introduction of 
HDMF in the TWRC model, a novel scheduling scheme should be 
proposed to effectively use its potential. 

We analyze the average queue lengths at the end nodes and relay 
and their instability. Firstly, define the following transmission modes 
for the whole system: 

• Mode I: A and B transmit data packets to the relay simultane¬ 
ously at the data rate of TZi and TZ 2 using the HDMF protocol; 

• Mode II: Relay broadcasts packets to A and B at the data rate 
of TZi if both of its queues have data; 

• Mode III: Relay only transmits packets to B using the data rate 
7 ^. 4 ; 

• Mode IV: Relay only transmits packets to A using the data rate 
7^5; 

A scheduling scheme is then proposed where the relay is set to be 
the administrative node. The scheduling details are given below: 

1) The Relay chooses a transmission mode from the above four 
modes with probability /i, (i = {1, 2, 3,4}). 

2) If Mode I is selected, both A and B will be instructed by Relay 
to transmit at the beginning of the time slot. If either of their 
buffers are not empty, e.g., min{(3a(f),Q(i(f)} > 0, the node 
with data in its buffer will transmit packets to the relay during 
this time slot; otherwise, they remain silent. 

3) If Mode II is selected. Relay will broadcast to A and B if both 
of their queues are not empty i.e., Qra and Qrb are greater 
than 0. In the case of one empty queue. Relay will broadcast 
the data from this queue. If both of them are empty. Relay will 
remain silent. 

4) The Relay will transmit to B if Mode III is selected and Qrb 
is not empty. If Mode IV is selected. Relay will transmit to A 
if Qra is not empty. If the queue for data to be sent to A or 
B at Relay is empty when the corresponding mode is selected. 
Relay will remain silent during this time slot. 


A. Queue Length Analysis 

Without loss of generality, we consider the route from A to B via 
Relay. The backward route from B to A through Relay is similar to 
the forward one and will be omitted here. As described before, the 
queue length of A is denoted by Qa and the queue length of Relay 
for the packets to be transmitted to B is denoted by Qrb- The state 
transition relation of {Qa, Qrb} can be modelled by a Markov chain 
using one-step state transition probability q{m,k},{i,j} which denotes 
the probability of the event that the state changes from Qa{t) = m, 
Qrb{t) = fe to Qa{t -f 1) = L Qrb(t -|- 1) = j. The probability that 
i packets arrive at one source node during the current time slot can 
be modelled by a Poisson distribution as follows 


where A is the data arrival frequency and T is the length of one time 
slot. An analysis method similar to El is introduced here. 

Under the condition of fading TWRC, channel quality can affect 
the data rate and queue length. The data rate which can be supported 
by the channel from A to Relay is denoted as n and the probability 
to support such rate is given by Similarly, the probability for 

Relay-A channel supporting a data rate of n is denoted as and 

the probability for Relay-B channel supporting n is qn{t). 

Under Mode I, with the application of HDMF protocol. Relay has 
two possible inputs: packets from differential DMF or direct DMF. 
Based on the analysis in Section EYl we can use Pabr to denote 
the probability with which the packets arriving from A (and B) by 
differential DMF, par to denote the probability of packets arriving 
from A by direct DMF and pbr to denote the probability of packets 
arriving from B by direct DMF. It is easy to know that pabr +Pabr -f 
Pbr = 1. The state transition probability is then given by 

I 

q{r. - 


(51) 





ai 

p 

tl 


m ^iT-iPahr + Par) 

0 < m < n, j 

= k + m 

ai-m 2^n=m ^ripbr 

0 < m < n, j 

= k 

ai — m-\-nan{jpabr H” Par) 

m > n > 0, j 

= k + n 

ai-mO-riPbr 

m > n > 0, j 

= k 


Notice that Mode I is selected with the probability of fi at the 
beginning of every time slot. The first item in | I5H denotes the state 
transition probability when the queue at A is empty: Qait) = m = 0. 
Therefore Qrb{t + 1) = Qrb{t) (which means j = k). The 
corresponding probability is m. 

However, if the queue at A is not empty, the state transition at Relay 
will depend on the HDMF protocol. If the channel can support a data 
rate n > m, the probability of receiving and storing m packets in the 
queue Qrb{t+1) (and thus causing the corresponding state transition) 
is aiEr=mCn(R abr+Par), because both differential DMF and direct 
DMF can generate the m packets. As a result, the queue length 
of Relay becomes j = k + m. The probability for the relay state 
remaining k is ai-m CnPbr, because Relay only receives data 

from B and these data would contribute to none of the two queues 
under concern. 

On the other hand, if n < m, only n packets will be received 
by Relay through either direct or differential DMF. Therefore, the 
probability of Relay’s state becoming k + n is ai t^a{Pabr + 

Par)- Another case is direct DMF for data from B, which would 
not contribute to the state transition. So j = k with a probability of 
tti — mCnPhr - 

Under Mode II - Relay broadcasts data to A and B, the state 
transition probability is given as follows 


II 

q{m.,k},{i,i] 


O-i — m ^2n — k 7 — 0 

ai-mV-n j = k — n,0 < n < k 


(52) 


If the channel can support a data rate of n which is greater than 
the number of packets within Relay’s buffer k, the whole buffer will 
be transmitted to the destination and thus emptied, e.g. j = 0. The 
probability for such state transition is ai-m E^fe ’’’n- Otherwise, the 
new queue length of Relay will become k — n, and the state transition 
probability is ai-mr-n- 

Under Mode III, Relay transmits packets to B and the state 
transition probability is as follows. 


III 

q{m.,k],{i,i] 


ai-m J2'^=k la i = 0, n > fc 
ai-mq-a j = k — n,G < n < k 


(53) 


If the data rate n equals or is greater than the data in the queue 
at Relay, j will be zero after transmission since all of its data has 
been transmitted. Otherwise, the queue length will be fc — n. The 
corresponding transition probabilities are given in the equation above. 


ai = 
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If Mode IV is selected, Relay will only transmit data to A. The 
state of our concerned queues will remain the same. Thus the state 
transition probability is as follows 

CLi — mj j k. (54) 

Notice that Relay selects the current transmission mode from the 
four candidates based on the probability of fi {i = By 

combining the one-step state transition probabilities and the mode 
choice probabilities, we can obtain the final transition probability, 
summarized in OS. The stationary state of the queue lengths of 
interest can be obtained given the one-step transition probabilities. 

In order to know the stationary-state queue length of Qa at A and 
Qrb at Relay, we need to work out the stationary-state distribution 
of the previous mentioned Markov chain with the transition proba¬ 
bilities summarized in l l52t . Firstly we rewrite the one-step transition 
probabilities using the matrix form. By denoting the buffer size for 
the queue at Relay for B as a big number N, the (x, y)th element of 
the transition matrix, e.g. is given by q{mN+k,iN+j} which 

is equivalent to q{m,k},{i,j}- The stationary-state distribution of the 
queue lengths can be obtained as follows 

n = l-(I-P + U)-\ (56) 

where 1 denotes a row vector with elements of ones, U is an all¬ 
unity matrix and I is the identity matrix 1261 15.107]. IT denotes the 
stationary probability of the two queue lengths {Qa{t),Qrb{t)} and 
can be used to calculate the average queue length of A, as follows 

OO OO 

Qa = 'y .'y , iT^in+j ■, (57) 

i=0 j=0 

and the average queue length of Relay for B, as follows 

OO OO 

Qrb ~ ^ ^ ^ j'^iri+j ■ (58) 

i — Q 0 

VI. Simulation 

A. Experimental Conditions 

We investigate the performance of the proposed HDMF protocol 
and the transmission scheduling scheme in this section. The chan¬ 
nels are modelled as block Rayleigh distributions and the noise is 
modelled as additive white Gaussian distributions for most of the 
cases, except the second experiment in Section IVl-BI where compare 
the protocols under Gaussian channel settings. We assume reciprocal 
channel settings, e.g. Har = Era and Err = Err- The modulation 
scheme used in the nodes is QPSK and each packet has 128 symbols. 
At every experiment, a total of 10® packets are transmitted from one 
source node to the other via Relay. 

B. Performance under Relative Channel Strengths 

The first experiment tests the Packet Error Rate (PER) perfor¬ 
mance of the HDMF protocol, the PNC protocol, the classic DMF 
protocol and the ANC protocol in Rayleigh fading scenarios. The 
average channel gains are adjusted using log(E{|/i_Bfl|}/E{|/iAH|)}. 
We set the range of this value from -0.7 to 0.7 and 0 means 
E{|/isij|} = E{|/iAi{|}. Figl^ shows that ANC and PNC have the 
highest PERs, which means their performance is heavily affected by 
the channel realizations. They also have their own minimum PER 
when both the two channels have the same average channel gains, 
e.g., E{|/is_r|} = E{|/iAfl|}. Classical DMF is able to achieve lower 
PER than ANC/PNC in this case. The proposed HDMF protocol 
follows a similar tread as ANC/PNC but with lower PER than all 
of them. Such trend becomes even clearer at high SNR settings, 
as shown in Fig[7] where the SNR is 25dB. In this scenario, the 
performance of all four protocols is improved; the most impressive 



Fig. 6. PER performance comparison of the protocols under Rayleigh channel 
settings and Eb/Ng = 15dB. 
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Fig. 7. PER performance comparison of the protocols under Rayleigh channel 
settings and Eb/Ng = 25dB. 

change is from PNC as its PER is below DMF HDMF has the best 
performance among them. 

The second experiment investigates the performance of these 
protocols under the Gaussian channel settings where the fading is 
small and merges into noise. Since the fading effect is avoided in 
this experiment, the impact of noise towards these protocols can be 
clearly identified. The channel gains vary by following the relative 
strength of log(|/i_BH|/|/iAi?|). A similar range is set from -0.7 to 
0.7, where 0 means \hRR\ = |/iAfl|. Figs. and [7] show the 

experimental results in terms of PER. 

Erom Eig. we can see that under Gaussian channel settings 
and SNR per bit Eb/No = lOdB, HDME and PNC achieve similar 
performance at the middle of the relative channel strength range 
( [—0.4, 0.4]), with its lowest PER at the point where the two channels 
have the same strength. Such advantages decrease when the channel 
gains of the two channels (A-Relay and B-Relay) diverge and finally 
they are overtaken by the DME protocol. As expected, DME has poor 
performance when the two channels have similar strengths because 
one channel heavily interferes the demodulation of the other. ANC 
has the worst performance because it amplifies noise significantly, 
especially when SNR is low. We raise SNR from lOdB to 15dB and 
obtain Eig. Under this new setting, the performance of ANC is 
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ai 

/ = fc = m = 0 

ai-m(l — /r + /iCo) 

j = k = 0,m > 0 

ai-m{f2rn + fsqn) 

j = k — n,k>n>0 

ai-m(h + U 9^) 

j = 0, k > 0 

aiifi + f 2 ro + foqo + U) 

j = k > 0,m = 0 

ai-miflPhr + /2f0 + fsqO + fi) 

j = k > 0,m > 0 

ai — m-Vnen(^Pabr 4” Par)fl 

j = k + n,m>n>0 

t^ri{Pabr + Par)fl 

j = k + m,m > 0 


0 else 



Fig. 8. PER performance comparison of the protocols under Gaussian channel 
settings, Ef,/No = lOdB. 



Fig. 9. PER performance comparison of the protocols under Gaussian channel 
settings, Ef,/NQ = 15dB. 

improved remarkably and it has lower PER than DMF. The PERs of 
HDMF and PNC are also improved significantly and have similarly 
low value at most of the x-axis range (They are 0 and not shown in 
the figure). 

C. Average SER 

This section compares the SER performances of the different 
relaying protocols. The experiment conditions are similar to the 
previous experiments, except for the following settings: 1) the average 



Fig. 10. SER performance of the four protocols. 

fading channel gains are fixed to 1; 2) Eh/No varies from 5dB 
to 30dB; 3) both the theoretical analysis and simulation use BPSK 
modulation. 

FigUo] demonstrates the average SERs of the four protocols. All 
of them have decreasing SER when SNR increases. However, at the 
lower SNR region (Eh/No < lOdB), DMF has lower SER than the 
other protocols. With the increase of SNR, the SERs of HDMF and 
PNC drop much faster than DMF. After 25dB, even ANC has a 
lower SER than DMF. For the regions with SNR > lOdB, HDMF 
outperforms PNC, ANC and DMF. 

FigHUcompares the average SER performance between simulation 
and theoretical analysis m. From the figure, we can see that 
simulation result matches well with theoretical analysis, particularly 
at high SNR regions for the reason of applying high SNR assumption 
in the theoretical analysis. Such results confirm the performance of 
the proposed protocol. 

D. Queue Length Analysis 

This experiment studies the queue length of A and the queue length 
for data to B at Relay in terms of Qa and Qrb- We assume the 
mode selection frequency fi to be equal for all the four modes, e.g., 
f{i) = 1/4, i = 1, ...,4. The channels are modelled as independent 
Rayleigh fading channels, e.g., Cn, Vn and qn follow the Rayleigh 
probability distribution function. Similar to da, the packet arrival 
frequency for A is A = 0.5 frames/slot. The time slot spans 1ms, e.g., 
T = 1ms, and the SNR is 20dB. We introduce the same design logic 
as in oa to configure the system to accommodate a higher average 
arrival rate than the actual rate of A :—>■ A(1 + e) for the purpose of 
reducing memory usage and reasonable running time, where e has a 
small positive value. 
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Fig. 11. Average SER comparison. 



Fig. 12. Stationary-state queue length of A and Relay for B versus e. Rayleigh 
fading channels are used. The two protocols HDMF and DNC ai‘e compared 
at the packet aniving rate of 0.5 frame/slot. 

Fig.im shows the stationary-state queue lengths for the proposed 
HDMF and the referenced scheduling for Digital Network Coding 
(DNC) d). It can be seen from the figure that the proposed 
HDMF has comparable smaller queue lengths both at Relay and 
node A than that of DNC. The reason is that HDMF utilizes the 
two advantages; reduced time slots consumption and opportunistic 
DMF under different channel conditions. LLR controls the protocol 
switching probabilities (i.e., Pabr, Par and pbr) to ensure the packet 
detection and effective packet forwarding. Generally, the queue length 
of HDMF is about 0.5 frames less than the DNC protocol. With the 
increase of e, the queue lengths of both protocols reduce significantly. 
As e is introduced to adjust the system throughput and capacity, 
however, the decrease in queue length occurs with the cost of 
decreased system throughput since there will be more idle slots during 
which no data is sent out. 

VH. Discussion 

This section discusses some issues involving the application of 
the HDMF protocol. The first issue is about multi-hop networks with 
several intermediate nodes. The demodulation and implementation of 
HDMF should be similar to the system model in this paper. However, 
data processing on relays should be carefully handled in order to 


avoid the same packet being sent back to its origin, thus reducing 
the efficiency. The routing algorithms developed for multihop sensor 
networks, e.g., (27), can be applied to solve this problem. (The 
specific implementation details are omitted in this paper.) 

This paper introduces the HDMF protocol with the default setting 
that the packets from two source nodes have the same length. If 
they are of different lengths, HDMF can still be applied with small 
amendments. Because the detection at Relay is based on LLR values, 
we can add some extra 0 bits to the shorter packet in order to calculate 
the LLR values. A better enhancement is to use these bits for the 
purpose of error-correcting coding 1251 . 

ANC and PNC can perform well if the two channels of TWRC have 
similar gains. This places a high requirement on synchronization as 
practical wireless channels cannot always guarantee the simultaneous 
arrival of signals from two distributed sources, thus the synchroniza¬ 
tion of Relay to one source would mean discrimination against the 
other and cause significant detection errors. If such case does happen, 
the proposed HDMF would be reduced to direct DMF automatically 
and avoid the negative influence from the severely discriminated 
channel which would degrades PNC or ANC significantly, and is 
thus less vulnerable to synchronization errors. 

The proposed HDMF does not try to solve the fundamental prob¬ 
lem of the signal to interference ratio being too low, rather it avoids 
the problem by relaying the signals through differential DMF rather 
than detecting the information forcefully which might be degraded 
by interference and noise. In this case, differential DMF provides 
more information for destination to be able to complete the final 
symbol detection, e.g. destination can reduce its own contribution in 
the interference. 

The instantaneous end-to-end error rate can provide a more accu¬ 
rate measure for the relay to make decision. However, its cost includes 
overheads incurred by feedbacks from the nodes and increased com¬ 
putation complexity. The next stage of the research will investigate 
these factors and try to establish an efficient solution based on the 
current protocol framework and the instantaneous end-to-end error 
rate. 

VIII. Conclusion 

This paper proposes an HDMF protocol for the TWRC model 
where existing protocols are troubled by high error rate or difficulty 
in synchronization. We study the components of HDMF and the 
fundamentals of direct DMF, differential DMF and the key detection 
criterion. We further analyze the queue length of the model with 
HDMF protocol. Comparison of its performance with the existing 
protocols and scheduling scheme indicates that the proposed HDMF 
has a lower average PER and smaller average queue length. 
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